Recovery of classical chaotic-like behaviour in a conservative quantum three body 

problem 



M.J. EveritlB 

Centre for Theoretical Physics, The British University in Egypt, El Sherouk 
City, Misr Ismalia Desert Road, Postal No. 11837, P.O. Box 43, Egypt. 

Recovering trajectories of quantum systems whose classical counterparts display chaotic behaviour 
has been a subject that has received a lot of interest over the last decade. However, most of these 
studies have focused on driven and dissipative systems. The relevance and impact of chaotic-like 
phenomena to quantum systems has been highlighted in recent studies which have shown that 
quantum chaos is significant in some aspects of quantum computation and information processing. 
In this paper we study a three body system comprising of identical particles arranged so that 
the system's classical trajectories exhibit Hamiltonian chaos. Here we show that it is possible to 
recover very nearly classical- like, conservative, chaotic trajectories from such a system through an 
unravelling of the master equation. Firstly, this is done through continuous measurement of the 
position of each system. Secondly, and perhaps somewhat surprisingly, we demonstrate that we still 
obtain a very good match between the classical and quantum dynamics by weakly measuring the 
position of only one of the oscillators. 
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I. INTRODUCTION 

Quantum mechanics is perhaps the most powerful and 
useful theory of physics to date. Indeed, with the pos- 
sible emergence of many new quantum technologies in 
areas such as computation, communication, cryptogra- 
phy and metrology this trend looks set to continue well 
into the future. With such strong interest in applica- 
tions of quantum mechanics comes a concomitant inter- 
est in the measurement process and the interaction be- 
tween ''classicaF and quantum systems. Indeed, as we 
wish to understand and apply quantum mechanics within 
the context of modern technology we will need to develop 
our understanding of what actually constitutes a classical 
device and how such objects interact with quantum sys- 
tems. However, the recovery of classical mechanics is not 
always as simple as a implied by a naivete interpretation 
of the correspondence principle. This essential require- 
ment of any physical theory can, for quantum mechanical 
systems, be stated as: 

'Tf a quantum system has a classical ana- 
logue, expectation values of operators behave, 
in the limit h ^ 0, like the corresponding 
classical quantities'^ [l| 

We observe that interpretation of this statement can be 
problematic if, for example, we consider quantum sys- 
tems that lack a specific dependence on Planks con- 
stant [2]. Further difficulties arise when attempting 
to recover the classical trajectories of classically non- 
linear and chaotic systems as the Schrodinger equation 
is strictly linear. 
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We note that these concerns are no longer just of in- 
terest to those studying either the measurement prob- 
lem or the correspondence principle and the emergence 
of the classical world. Indeed this area has a direct im- 
pact on quantum technologies. In order to fully leverage 
the power afforded by these emerging fields we must not 
only understand in depth the measurement process but 
also many body quantum systems coupled to real envi- 
ronments. This is highlighted by the recent observation 
of chaos in the spectrum of Shor's algorithm [3] as well 
as in other studies involving quantum information pro- 
cessing and quantum chaos 0, S @] . 

A solution to the correspondence problem for chaotic 
systems which has been employed with great success is 
found by utilising quantum trajectories methods 0,0,0, 
M, EH, E, E, H, IB, [H, El- Here, introduction of en- 
vironmental degrees of freedom and unravelling the mas- 
ter equation yield stochastic Schrodinger equations from 
which chaotic-like trajectories may be recovered. This 
process can be considered as comprising of several steps. 
Firstly, we make the quantum system of an open one. 
This is archived by coupling the quantum object to an 
environment which may take the form of a measurement 
device. Once the environment has been introduced we 
model the evolution of the system's density operator (in 
the presence of the environment) using a linear master 
equation. However, master equations are similar to the 
Langevin equation insofar as they only predict a set of 
probable outcomes over an ensemble of systems or ex- 
periments. Therefore, in order to get some idea of the 
possible behaviour of an individual experiment we next 
unravel the master equation. In essence, this process in- 
volves finding a stochastic differential equation for the 
system's state vector with the proviso that the dynamics 
given by the master equation are returned in the ensem- 
ble average over many solutions. There are an infinitely 
many ways to do this each representing a different phys- 
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ical process. In this work we employ the quantum state 
diffusion (QSD) unravehing which corresponds to a unit- 
efficiency heterodyne measurement (or ambi-quadrature 
homodyne detection) on the environmental degrees of 
freedom [16] (for an detailed introduction to this ap- 
proach please see [HI). Here the evolution of the state 
vector is given by the Ito increment equation [l3l. fl^ 
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where the Lindblad operators Li represent coupling to 
environmental degrees of freedom, dt is the time incre- 
ment, and (i^ are complex Weiner increments defined 
such that = ^ = and d^ = dt [H, 
Throughout this work a bar over a quantity denotes the 
average over stochastic processes whilst the notation (•) 
is used for quantum mechanical expectation values. The 
first term on the right hand side of this equation is the 
Schrodinger evolution of the system while the second 
(drift) and third (fluctuation) terms describe the deco- 
hering effects of the environment on the evolution of the 
systems state vector. 

However, to date the body of work which uses quan- 
tum trajectories to recover classically chaotic-like trajec- 
tories has focused on those systems that are dissipative. 
There are several notable exceptions that demonstrate 
that continuous measurement of both driven and un- 
driven conservative systems that can recover classical like 
behaviour [l^, [2O, [U, [22[ . However, these works consider 
systems with only a single degree of freedom. Recently 
we became interested in whether it was possible, using a 
similar analysis, to recover chaotic trajectories of classi- 
cal, multi-component, systems undergoing Hamiltonian 
chaos. Initially we wished to consider the traditional 
three body problem of classical mechanics for particles 
with similar masses [30]. This problem, although his- 
torically very significant, is non- trivial to solve. Conse- 
quently, in this work we consider a somewhat simplified 
system comprising of three coupled one-dimensional an- 
harmonic oscillators. 



II. BACKGROUND 

The Hamiltonian for our chosen three body system, 
comprising of one-dimensional anharmonic oscillators 
with a quart ic potential and unit mass, is given by 
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The classical dynamics associated with this Hamiltonian 
can be chaotic and are known to have positive Lyapunov 



exponents [23|. When we consider classical mechanics 
Qi and Pi are taken to represent the classical values of 
position and momentum. However, when we consider the 
quantum mechanics they are replaced by their operator 
counterparts. As we shall always be clear as to which 
description we are considering at any one time this does 
not lead to any ambiguity. 

We have already stated one expression of the corre- 
spondence principal in quantum mechanics. An alterna- 
tive definition, which we find preferable, is to 

^^consider h fixed (it is) and scale the Hamilto- 
nian so that the relative motion of the expec- 
tation values of the observable become large 
when compared with the minimum area {h/2) 
in the phase space". 

In either case this is the role of the (3 term in the Hamil- 
tonian, i.e. h ^ ph so that the smaller P the larger the 
dynamics when compared to a plank cell. 



III. RESULTS 

From the Hamiltonian ([2]) we find the three classical 
equations of motion are 




0, where z, j = 1, 2, 3. (3) 



When we solve these coupled equations of motion with 
the initial conditions qi = — 0.2//3, q2 = 0.05//3, qs = 
0.15 /P and pi = for all i we find that the dynamics 
are chaotic. We show the phase portrait for the solutions 
to these equations in figure [11(a) where, without loss of 
generality, we have set P = 1. We note as an interesting 
aside, one feature of this system is that if qi = pi = for 
any i then qi = pi = always. 

We now proceed to discuss the quantum mechanical 
description of these coupled oscillators. Unlike the classi- 
cal equations of motion the Schrodinger equation for this, 
or any other, system is strictly linear. We find that so- 
lutions to the Schrodinger equation for this system, even 
for moderate values of /3, delocalise so rapidly that ob- 
taining accurate solutions is not possible for us. This 
does, however, reinforce the lack of correspondence be- 
tween classical dynamics and Schrodinger evolution for 
this system. 

Following past work 0, i, i, [O, E, E B E 
\Tdi . fTrt on recovering classically chaotic like orbits from 
a system's quantum counterpart we solve the unravelling 
of the master equation ([1]) with Hamiltonian ([2]). For 
this example there are three points of note with regard 
to possible choices of the environmental degrees of free- 
dom. Firstly, coupling to an environment helps localise 
the system's state vector and hence produce a well de- 
fined, classical-like, trajectory. Secondly, as the classi- 
cal system is Hamiltonian and therefore conservative, we 
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FIG. 1: (colour on-line) An example chaotic trajectory (a) 
solutions to the classical equations of motion ^ for the ini- 
tial conditions qi = —0.2, q2 — 0.05, ^3 = 0.15 for our chosen 
system and (b) scaled quantum expectation values f3 (qi) and 
P (pi) versus time for an unravelling of the master equation 
with initial state L>i(-0.2//3)L>2(0.05//3)L>3(-0.15//3) |000) 
and f3 = 1/2000. Trajectories for oscillator one are shown 
in magenta (medium grey), for two in blue (dark grey) and 
three in green (light grey). Note, all quantities are dimen- 
sionless. 
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FIG. 2: (colour on-line) A comparison between the classi- 
cal position qi (dashed grey) and f3 (qi) for three different 
couplings (k) to the environment. Note, all quantities are 
dimensionless. 



must chose the environment of each oscillator so that en- 
ergy exchange is minimised between any part of the sys- 
tem and the environmental degrees of freedom. Thirdly, 
we should specify a physically reasonable environment. 

In this work we have chosen initially one of the most 
obvious candidates for the environment which satisfies all 



these conditions. Explicitly we have set each Lindblad 
Li = /^^i, i = 1,2,3 corresponding to the continuous 
measurement of position. This unravelling also corre- 
sponds to that of the Master equation for a weakly cou- 
pled, high temperature, thermal environment [l^. Here 
hi represents the magnitude of the coupling between each 
component of the system and its respective environment. 
In this work we use several values of hi. In figure WO^) 
we use an intermediate coupHng {tz = 0.1) whilst in 
figures El [3l [H and [5] we also present results for weak 
{hi = 0.01) and strong {hi = 0.5) couplings. 

As our initial state, and for the best possible match 
with the classical initial conditions, we chose a tensor 
product of coherent states for which the quantum ex- 
pectation values in position and momentum are centred 
in q — p phase plane at qi = — 0.2//3, = 0.05//3, = 
0.15//3 and Pi = where i = 1, 2, 3. Alternatively, we can 
express this initial condition explicitly as translated vac- 
uum states by L>i(-0.2/^)L>2(0.05//3)L>3(-0.15//?) |000) 
where Di{.) is the displacement operator in position for 
each component of the system. Here we have chosen 
(3 = 1/2000 as this is the smallest value for which we can 
solve ^ both accurately and within a reasonable time 
frame. In order to help the reader quantify the time 
scale over which our results are presented we note that 
the log time associated with our chosen value of beta is 
log(l//3) ^ 3.3 [2i,[25[. This is much shorter than the pe- 
riod over which we present the evolution of the system's 
trajectories. 

In figure [Ifb) we show the dynamics of the quantum 
expectation values of the position and momentum op- 
erators for each oscillator. These have been scaled by 
a factor of /3 so that they may be compared with fig- 
ure [TJa). Here we see very good agreement initially, and 
similar characteristics throughout, the displayed dynam- 
ics. Indeed the trajectories are similar enough that it is 
impossible to determine from the graph alone which plot 
shows the classical and which the quantum evolution. We 
note that these curves begin to differ after a short period 
of time. However, this is not unexpected as the system 
we are analysing is chaotic. In order to make the rea- 
sonable comparison of these results readily available we 
also include a graph of the evolution of qi and {qi) as a 
function of time in Fig. [2] for three different couplings to 
the environment. It is apparent that there is a very good 
match between the quantum expectation values and the 
classical trajectory for hi = 0.01. 

We note that simply by including an environment our 
system no longer undergoes Hamiltonian evolution. In 
other words, in order to be able to recover classical like 
trajectories of quantum systems whose classical counter- 
parts exhibit Hamiltonian chaos we include environmen- 
tal degrees of freedom that imply non-Hamiltonian evolu- 
tion of the quantum system. However, using a sufficiently 
low coupling strength to the environment results in a con- 
comitant reduction both in energy exchange between the 
system and it's environment and the localisation of the 
state vector. We now verify that the solutions to equa- 
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FIG. 3: (colour on-line) Total system energy, with magni- 
fied section inset, computed by substituting (qi) and (pi) into 
the Hamiltonian (|2| for three different couplings (k) to the 
environment. Note, all quantities are dimensionless. 



tion ([T]) are to, good approximation, conservative. This 
does indeed appear to be the case for both intermediate 
and weak couphngs {k, = 0.1 and 0.01) but not for the 
stronger couphng {tz = 0.5). This can be seen in figured 
where we show the total energy found by substituting 
{qi) and (pi) for q^ and pi into the Hamiltonian ([2]), i.e. 



Energy : 
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E (4) 



where z,j = 1,2,3. We note that we do not compute 
{H) as we wish to compare directly with the equivalent 
classical calculation. 

Next we verify localisation of the state vector by com- 
puting the uncertainty in position and momentum for 
the first oscillator for three different couplings to the 

environment. Because both Ag^ = (qf) — {qi)'^ and 



^Pi — \l {Pi) ~ (Pi) between components behave in a 
similar fashion, we do not show results for the other two 
oscillators here. As is evident from figure [H the interac- 
tion with the environment causes system's state vector 
to localise within each of the component spaces. It is 
also apparent from this figure that the level of localisa- 
tion is dependent of the coupling between each of the 
system's components and their respective environments. 
We also note that as the system evolves it's states become 
squeezed in each of the momentum variables. Unlike the 
results presented in figure [TJb) we have not scaled these 
uncertainty values by = 1/2000. Consequently, for di- 
rect comparison with the first two figures the results pre- 
sented in figure m should be divided by 2000. Hence, the 
uncertainty in either position or momentum can be seen 
to be quite negligible when compared with the trajectory 
of quantum expectation values plotted in figure [Hb). 

We can extract further information on the dynamics of 
this system simply by borrowing a technique from quan- 
tum optics. Namely through analysing the photon statis- 
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FIG. 4: (colour on-line) Uncertainty in position (light 
grey /magenta) and momentum (dark grey/blue) as a function 
of time for the first component for three different couplings 
(k) to the environment. Values beneath the dashed (green) 
line indicate squeezing. Note, all quantities are dimensionless. 



tics (bunching of photons) described by the second order 
correlations [261, \2l\ 



(2) _ (nf ) - {rii) 



1,2,3 



(5) 



Where n is the number operator. Values of g^'^^ greater 
than 1 indicate photon bunching where photons arrive 
in groups while values of g^'^^ smaller than one indicate 
antibunching, a purely quantum mechanical phenomena 
representing precisely regular arrival of photons. How- 
ever g^'^^ = 1 corresponds to Poissonian statistics and 
which is what we would expect from the state of our sys- 
tem should it be undergoing a classical like evolution. As 
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FIG. 5: gf^ (i = 1) coefficient, with magnified section inset, 
as a function of time for tfiree different couplings (k) to tfie 
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environment. Classical- like motion yields g, 
implies a Poissonian statistics for the state of the system. 
Note, all quantities are dimensionless. 
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FIG. 6: (colour on-line) An example chaotic-like trajectory 
for the scaled quantum expectation values f3 (qi) and f3 (pi), 
for weak measurement of the first oscillator only. Here Li = 
Kqi,L2 = Ls = for = 0.01 and p 1/2000. Quantum 
trajectories for oscillator one are shown in magenta (medium 
grey), for two in blue (dark grey) and three in green (light 
grey). The corresponding classical dynamics are shown with 
a dashed grey line. Note, all quantities are dimensionless. 



we can see from figure [5] this is indeed the case (where, 
again, we have only shown data for the first component). 

In the discussion above we have considered what hap- 
pens when we perform simultaneous, continuous, mea- 
surements of the position of each of the components 
of the system. As the weak measurement limit is ap- 
proached we find good agreement between classical and 
quantum dynamics and the system becomes, to good ap- 
proximation, conservative. We now demonstrate that 
there exists a weaker condition under which we can pro- 
duce the same outcome. That is, we show that weak 
measurement of only one of the position variables is suf- 
ficient produce, to very good approximation, classical- 
like trajectories. As no one oscillator has a privileged 
status over the others we set, without loss of gener- 
ality, Li — nqi^L2 = Ls = with = 0.01 and 
f3 = 1/2000. Again we solve the unravelling of the 
master equation ^ with Hamiltonian ([2]) and the ini- 
tial state Di{-0.2/ p)D2{0m/ p)D3{-0.15/ p) |000). As 
we can see from figure [6] even under these conditions we 
recover quantum trajectories whose expectation values 
match very well indeed with those of the equivalent clas- 
sical system. 



IV. CONCLUSION 

For any given /3 and an initially localised state there 
will be some agreement between the dynamics of the clas- 
sical system and the evolution of the quantum expec- 
tation values of the corresponding quantum operators. 
However, after a short period of time the quantum state 
vector will begin to delocalise and differences between 
the predictions of each theory become apparent. We 
have demonstrated, by localising the state vector through 
measurement of the position of one or all components of 
the system, that near classical like dynamics can be re- 



covered through an unravelling the master equation. For 
this work we have chosen quantum state diffusion. How- 
ever, it is likely that any other unravelling will produce 
similar results. Such detailed analysis is beyond the scope 
of this work and would belong in a more in depth study. 
We also note, as a subject for future study, that it would 
be interesting to determine the conditions under which 
measurement of a subset of the degrees of freedom of 
an A^-body system results in the localisation of the state 
vector. 

Finally we would like to observe that following [2^ it 
would be interesting to characterise the entanglement be- 
tween the components of this system. As it may well be 
the case that for this example, as well as the one studied 
in [28] , that the entanglement does not necessarily vanish 
in the classical limit. Unfortunately current restrictions 
on computational power prevent us from conducting this 
study at the present time. However, from the last re- 
sult presented here we intuitively feel that there must 
persistently exist at least a small degree of entanglement 
between the first component and each of the other two. 
In order to justify this statement we propose the follow- 
ing argument. First, consider the extension to the tensor 
product space of the Lindblad operator, explicitly this is 
tzqi I2 I3. Now, by examining equation ^ we see 
that if the state of the system was separable, the last two 
terms (those responsible for the localisation of |?/^)) would 
not affect the components of the state vector for the sec- 
ond and third degrees of freedom. As introduction of 
this environmental degree of freedom localises the state 
vector and result in the recovery of a classical like tra- 
jectory, we therefore propose that the state vector must 
posses some non-zero entanglement. 
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